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Self-accelerating solutions in massive gravity provide explicit, calculable examples that exhibit the 
general interplay between superluminality, the well-posedness of the Cauchy problem, and strong 
coupling. For three particular classes of vacuum solutions, one of which is new to this work, we 
construct the conformal diagram for the characteristic surfaces on which isotropic stress-energy per¬ 
turbations propagate. With one exception, all solutions necessarily possess spacelike characteristics, 
indicating perturbative superluminality. Foliating the spacetime with these surfaces gives a patho¬ 
logical frame where kinetic terms of the perturbations vanish, confusing the Hamiltonian counting 
of degrees of freedom. This frame dependence distinguishes the vanishing of kinetic terms from 
strong coupling of perturbations or an ill-posed Cauchy problem. We give examples where spacelike 
characteristics do and do not originate from a point where perturbation theory breaks down and 
where spacelike surfaces do or do not intersect all characteristics in the past light cone of a given 
observer. The global structure of spacetime also reveals issues that are unique to theories with two 
metrics: in all three classes of solutions, the Minkowski fiducial space fails to cover the entire de 
Sitter spacetime allowing worldlines of observers to end in finite proper time at determinant singu¬ 
larities. Characteristics run tangent to these surfaces requiring ad hoc rules to establish continuity 
across singularities. 


I. INTRODUCTION 

Building upon earlier work [1-3] , de Rham, Gabadadze 
and Tolley (dRGT) [4] first constructed a consistent in¬ 
teracting theory of a massive spin-2 field using an aux¬ 
iliary flat fiducial metric. Though the study of this the¬ 
ory and other nonlinear massive gravity variants has now 
reached a somewhat mature stage [5] , there remain many 
outstanding questions to be addressed and understood. 

Foremost amongst these is the status of cosmologi¬ 
cal solutions within the theory. The dRGT theory does 
not admit flat or closed cosmological solutions where the 
spacetime and fiducial metrics are simultaneously homo¬ 
geneous and isotropic [6, 7]. Though this complicates 
matters computationally, it does not automatically de¬ 
stroy the phenomenological viability of cosmological so¬ 
lutions. The spacetime upon which matter fields prop¬ 
agate can be homogeneous and isotropic for any choice 
of spatial curvature [8-10]. Indeed, many explicit exam¬ 
ples of cosmological solutions to the theory have been 
found where accelerated expansion occurs without a true 
cosmological constant [6-16]. 

Though these solutions are satisfactory as far as back¬ 
ground evolution is concerned, the theory of fluctuations 
about these solutions is somewhat paradoxical. Analy¬ 
ses of perturbations about these self-accelerating back¬ 
grounds [17, 18] reveal pathologies such as unbounded¬ 
ness of the Hamiltonian for isotropic perturbations, or 
strong-coupling of degrees of freedom in specific back¬ 
grounds, but often these appear in quite a subtle way 
[19]. For example it has been claimed that the number 
of degrees of freedom identified by a Hamiltonian analy¬ 
sis depends on a choice of coordinates [18]. Therefore, a 


worthwhile aim is to understand more fully the relation¬ 
ship between these potentially problematic features. 

More generally, this broad class of self-accelerating so¬ 
lutions provides an interesting playground to explore fun¬ 
damental field-theoretical questions which remain unre¬ 
solved in massive gravity and beyond. In particular, 
there has been much recent interest in the interplay be¬ 
tween superluminality in field theories and strong cou¬ 
pling phenomena, as well as what if anything the presence 
of these features implies about the well-posedness of the 
Cauchy problem (for various perspectives, see [20-39]). 
The self-accelerating sector of dRGT allows us an op¬ 
portunity to explore these issues with explicit solutions: 
the theory of isotropic perturbations about these solu¬ 
tions turns out to be rather simple [17]. In particular, it 
is possible to solve exactly for the characteristic hyper¬ 
surfaces upon which graviton stress energy fluctuations 
propagate. This makes it rather easy to study pertur¬ 
bative aspects of the causal structure of self-accelerating 
backgrounds. Therefore, another worthwhile avenue to 
explore is to what extent we can abstract general lessons 
about superluminality, strong coupling and Cauchy sur¬ 
faces from these particular examples. 

In this paper, we pursue both of these goals in tandem. 
After reviewing the structure of the dRGT theory—and 
in particular the construction of self-accelerating solu¬ 
tions of Ref. [8]— we focus on three particular families of 
vacuum solutions, one of which is new to this work. We 
explicitly solve for the characteristics for isotropic per¬ 
turbations and examine their causal structure via con¬ 
formal diagrams. In this sense, our study is related to 
the analysis of Ref. [35-40] but our exact solutions are 
general and do not require discontinuous or “shock” con- 
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ditions which automatically entail strong coupling. We 
find many interesting phenomena: in particular, the so¬ 
lutions we consider exhibit superluminal propagation of 
fluctuations by necessity except in one unique case, only 
some of which provide a well-posed Cauchy problem and 
originate from singular conditions where strong-coupling 
might reside. 

Perturbative superluminality on specific backgrounds 
does not necessarily indicate superluminality in the full 
theory (see [21-23]), which would present problems for a 
local and Lorentz-invariant UV completion of the theory 
itself [ 20 ], nor does it imply acausal structures such as 
closed timelike curves [41, 42]. However these examples 
do highlight the fact that the highly related notions of 
superluminality, strong coupling and well-posed Cauchy 
problem are conceptually distinct. 

Our examples also highlight an issue that is unique to 
a theory with two metrics. In all three families of vac¬ 
uum solutions, the Minkowski fiducial space fails to cover 
the entire spacetime. The point at which the Minkowski 
chart, or unitary gauge, ends is called a determinant sin¬ 
gularity; it is diffeomorphism-invariant and hence phys¬ 
ical [43, 44]. In these examples, worldlines of some ob¬ 
servers can intersect the singularity in finite proper time 
and for others the singularity lies in their past light cone 
requiring ad hoc rules for continuing worldlines or bound¬ 
ary conditions. 

The paper is organized as follows. In §11, we review the 
construction of bi-isotropic self-accelerating background 
solutions [ 8 ] and isotropic perturbations around them 
[17]. We then specialize to the vacuum case in §111 and 
discuss the global structure of the spacetime background 
and fiducial metric for three families of solutions, one of 
which is new to this paper. In §IV, we employ a char¬ 
acteristic analysis of perturbations around these back¬ 
grounds to expose the relationship between superluminal¬ 
ity, strong coupling and the well-posedness of the Cauchy 
problem. We discuss the implications of these results in 

§v. 

II. SELF-ACCELERATION IN MASSIVE 
GRAVITY 

In this section we briefly review the properties of the 
dRGT theory of massive gravity, the construction of self- 
accelerating background solutions, and spherically sym¬ 
metric perturbations around them. 

A. Fiducial Metric 

The dRGT [4] nonlinear theory of a massive spin-2 
field is given by the following action which propagates 
only the expected 5 polarizations of a massive graviton: 

S = ~ m2 y Fk ’ w 


where 1 = ( 87 rG ) _1 is the reduced Planck mass and 
the Ff~ terms are characteristic polynomials of the matrix 
7 . These can be written explicitly as 

Foh) = 1 , 

^1(7) = W, 

^2(7) = [7]M7 2 ], ( 2 ) 

^3(7) = [7] 3 -3[ 7 ][7 2 ]+2[7 3 ], 

^ 4 ( 7 ) = [ 7] 4 - 6[ 7 ] 2 [7 2 ] + 3 [ 7 2 ] 2 + 8[ 7 ][7 3 ] - 6 [ 7 4 ], 

where [] denotes the trace of the enclosed matrix. The 
matrix 7 is the square root of the product of the inverse 
spacetime metric g _1 and a flat fiducial metric £ 

7 ^ 7 “ = ■ ( 3 ) 

The parameters of the theory defined by Eq. (1) are 
m —the graviton mass—and the /3k- These parameters 
are not all independent, but rather they depend on two 
fundamental independent parameters { 0 : 3 , 04 } through 

A) = —12(1 + 2a 3 + 2aq), 

Al = 6(1 + 3a 3 + 4aq), 

ft =-2(1 + 6 a 3 + 12- a 4 ), (4) 

As = 6 (a 3 +4a 4 ), 

At = —24a 4 . 

The chart of the spacetime metric for which the fiducial 
metric is represented by the Minkowski metric = r)^ 
is called unitary gauge and there the spacetime met¬ 
ric possesses additional degrees of freedom compared to 
Einstein gravity. Diffeomorphism invariance can be re¬ 
stored by employing the Stuckelberg trick to introduce 
four scalar functions, <j> a , to represent the flat fiducial 
metric in arbitrary coordinates 

= Va-bd^d^. ( 5 ) 

The <fi a are equal to the unitary gauge coordinates x£ and 
absorb the extra polarization states in arbitrary gauges 
where diffeomorphism invariance is used to eliminate or 
constrain the metric terms. The matrix d^if> a then repre¬ 
sents the Jacobian of the coordinate transform between 
a general set of coordinates x M and x£. 

Minkowski space may not cover, or equivalently uni¬ 
tary gauge may not chart, the entire spacetime; this 
situation is signaled by a non-invertible Jacobian trans¬ 
form, which we call a determinant singularity [43]. Unlike 
a pure coordinate singularity, a determinant singularity 
does not depend on the chart of the spacetime. This is 
because the ratio of determinants of the two metrics 

det(g^ 1 S) = det(g” 1 r 7 ) = -det(g^ 1 ) (6) 

is a spacetime scalar. A coordinate singularity in the 
unitary chart of the spacetime metric g M therefore be¬ 
comes a coordinate-invariant determinant singularity. In 
contrast to a curvature singularity, the two metrics need 
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not individually have any diffeomorphism invariant sin¬ 
gularities, despite the fact that combined they display a 
determinant singularity. 

Physically, the presence of two metrics means that 
worldlines of observers may end after a finite proper time 
in spacetime has elapsed, since an infinite interval can 
elapse as measured by the fiducial metric. This geodesic 
incompleteness is another indicator that we should take 
determinant singularities seriously. It is possible to con¬ 
tinue worldlines past these singularities in the space- 
time, but this requires multiple copies of Minkowski space 
and a rule for continuing the chart—or equivalently the 
Stiickelberg fields—that is not directly imposed by the 
action. 


B. Background Equations of Motion 

For any isotropic spacetime metric, including cosmo¬ 
logical solutions with arbitrary matter content, there are 
solutions to the dRGT equations of motion where the 
stress-energy associated with the graviton potential in 
Eq. (1) behaves as a cosmological constant. The con¬ 
struction of Ref. [8] , which we now review, is in fact ex¬ 
tensible beyond dRGT to cases with non-flat bi-isotropic 
and dynamical metrics [44, 45]. 

Any metric with spatial slices invariant under SO(3) 
rotations can be written in isotropic coordinates, in which 
the line element takes the form 

ds 2 = -& 2 (f,r)dt 2 + a 2 (f,r)(dr 2 + r 2 dU 2 ), (7) 

where dfl^ is the h ne element on a 2-sphere. If the fiducial 
metric is rotationally invariant in the same coordinate 
system (be., the metrics are bi-isotropic), the Stiickelberg 
fields take the following form, 

4 >° = 

P= 9 (t,r )-, (8) 

r 

and are completely specified by the two functions / and 
g. Note that / = t u is the unitary gauge time and g = r u 
is the unitary gauge radius that describe the fiducial flat 
line element 

ds| = E^daJMaA = —d/ 2 + d g 2 + g 2 dill- (9) 

Unitary gauge uniquely specifies the coordinates, 
whereas the isotropic condition does not, leading to mul¬ 
tiple paths to finding the same solution and superficially 
different descriptions of their dynamics. Once solutions 
for / and g are found using any isotropic construction 
they may be re-expressed in an alternate choice of coor¬ 
dinates since both are spacetime scalars. 

Upon inserting the ansatze (7) and (8) into the ac¬ 
tion (1), we find that the equation of motion for the spa¬ 
tial Stiickelberg fields is satisfied by 

g(t,r) = x 0 a(t,r)r, (10) 


where the constant Xq solves the polynomial equation 
Pi(a;o) = 0 with 

Pi{x) = 2(3—2x)+6(®—l)(a;—3)a 3 +24(a;— l) 2 a 4 . (11) 


On this solution, the effective stress tensor due to the 
presence of the non-derivative graviton interactions takes 
the form of an effective cosmological constant 


I )iv = —(12) 


where 

A e ff = ^m 2 P 0 (x 0 ), (13) 

and the polynomial Pq(x) is given by 

Pq(x ) = —12 — 2x(x — 6) — 12(cc — l)(x — 2 )a 3 

- 24(a; - l) 2 a 4 . (14) 

Unitary time, f(t,r), satisfies the equation 
/— W 

v A' =-h x 0 , (15) 

X 0 

where 


A = 


- (f 




(J , 9 y (9 , J \ 

\b + »- a ) -\b + »Y) 

{jsf - gf), 


(16) 


with branches due to the matrix square root 7 defined in 
Eq. (3) allowing g = ±1. Here and throughout, over¬ 
dots denote derivatives with respect to t and primes 
denote derivatives with respect to r. Note that W is 
proportional to the determinant of the Jacobian trans¬ 
form from unitary gauge to isotropic coordinates. When 
W = ± 00 , 0 or is undefined because either / or g are not 
continuously differentiable, the Jacobian transform is not 
invertible. We call all of these cases a determinant singu¬ 
larity. Analytic continuation is sometimes possible, espe¬ 
cially in the later two cases, but requires a second fiducial 
metric and solution with its own choice of branch [43]. 
We pick g = 1 for the examples in the following sections. 

Using Eq. (10) in Eq. (15), we see that the latter equa¬ 
tion can be cast as 

b 2 f' 2 + 2ar(a! f 2 - aff) + r 2 (a'f - af) 2 

= x 2 0 ( a ,2 b 2 r 2 + 2a ab 2 r - dVr 2 ) . (17) 

This nonlinear partial differential equation has an infi¬ 
nite number of distinct self-accelerating solutions, each 
of which possesses the same background spacetime met¬ 
ric and A eff . In order to specify a solution, one assumes 
a functional form for a and b consistent with the Ein¬ 
stein equations sourced by A e fj and the matter content 
and then solves Eq. (17) to determine /. In §111, we will 
perform this procedure, choosing a and b so that there 
is no matter content and the background spacetime is de 
Sitter. 
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For solutions to the equation of motion (17) which sat¬ 
isfy the condition 

if = gg', (is) 

the fiducial metric Eq. (9) is also diagonal in the same 
(t, r ) coordinates that the spacetime metric is diagonal. 
We shall see that this is exactly the condition for which 
the kinetic term of isotropic perturbations vanishes [18]. 
However, we can analyze the dynamics of the same solu¬ 
tion in alternate isotropic coordinate systems which mix 
the temporal and radial coordinates where bi-diagonality 
and vanishing kinetic terms do not apply. Moreover, the 
diffeomorphism invariance of the Stiickelberg represen¬ 
tation means that we are not even limited to isotropic 
coordinate systems when analyzing the dynamics. 


C. Isotropic Perturbations 

Given the rotational invariance of the background 
spacetime metric and Stiickelberg fields, isotropic pertur¬ 
bations about these background solutions are straightfor¬ 
ward to analyze (c/. [19])- This is a reasonable starting 
point, as uncovering problems in this restricted class of 
perturbations would already indicate a pathology of the 
background solution. Note that the converse is not true: 
a background solution with healthy isotropic perturba¬ 
tions can still show pathologies in the anisotropic sector. 

In order to study fluctuations about the background 
solutions we are considering, we perturb both the metric 
variables and Stiickelberg fields about their background 
values; isotropic perturbations are specified by four func¬ 
tions of (t,r): 5a, 5b, 5f and 5g. Varying the quadratic 
action for these variables given in Ref. [17] yields an in¬ 
dependent equation of motion for a specific combination 

5Y(t,r) = 5g(t,r) — Xor5a(t,r). (19) 

This combination is special because it is precisely the 
variable which quantifies perturbations away from the ef¬ 
fective cosmological constant solution; perturbations for 
which <5T = 0 still satisfy (10) and therefore leave A e ff 
unchanged. The equation of motion for the variable <5T 
is [17] 


( 20 ) 

Equation (20) is first order in both time and space deriva¬ 
tives and does not depend on other perturbations. The 
5f equation of motion is also first order, but unlike the 
<5T equation it requires specification of <5T and other per¬ 
turbations. To linear order in the Stiickelberg fluctua¬ 
tions, the stress energy tensor of the perturbations de¬ 
pends only on AT and we therefore focus on its dynam¬ 
ics. The equation of motion (20) is equivalent to local 
conservation of stress energy [17]. 




, 2 2 
+ga r 


b ' a 
(ary' 


rJ 

= d r 

ST- 




abr ( g f' 


The coefficient of the temporal derivative term, 

A t = - fj,ar(ar)', (21) 

is of special interest to the time evolution of the perturba¬ 
tions. In a choice of isotropic coordinates where A t = 0, 
the kinetic term for 5T vanishes and hence the equation 
of motion (20) does not specify the evolution of perturba¬ 
tions off of a constant time surface. Combining Eq. (17) 
and Eq. (18), we see that A t = 0 whenever the spacetime 
and fiducial metrics are bi-diagonal and vice versa. 

When A t = 0, the energy density, momentum and 
pressure associated with any 5F fluctuation vanish to 
linear order [17]. Thus energy-momentum conservation, 
though enforced by the equation of motion, also does not 
yield equations that allow evolution of <5T off of constant 
time surfaces. Instead, the nonvanishing anisotropic 
stresses must obey a constraint equation. 

On the other hand, the fact that the equation of 
motion can be written in terms of different choices of 
isotropic time—or more generally different foliations of 
the spacetime—indicates that the vanishing of the ki¬ 
netic term is a coordinate dependent statement [18]. 1 In 
a choice where the kinetic term is present, the equation 
of motion does supply an evolution equation, or equiv¬ 
alently the constraint on stresses becomes a non-trivial 
equation for energy-momentum conservation. We shall 
see that this is a generic property for perturbations that 
propagate superluminally, i.e., on spacelike characteris¬ 
tics. 

Relatedly, Ref. [18] show through a Hamiltonian anal¬ 
ysis that the number of propagating (isotropic) degrees 
of freedom in a given isotropic frame is one if A t ^ 0 and 
zero if A t = 0. The number of degrees of freedom counted 
in this way, therefore depends on the coordinate system. 
For At = 0, the Hamiltonian is not associated with the 
time evolution of the system, in the sense that it defines 
evolution along a spatial slice, rather than transverse to 
it. As such, it is somewhat unclear how to interpret the 
Hamiltonian analysis in this case. Further, in the case 
At 7 ^ 0 , the canonical momentum for <5r appears lin¬ 
early in the Hamiltonian, and thus the Hamiltonian is 
unbounded from below. 


III. VACUUM SOLUTIONS 

In this section we specialize to background solutions 
without any matter content. In these solutions the only 
source of stress-energy determining the background is the 
effective cosmological constant, A e ff, and the spacetime 


1 Note that there are more prosaic examples of this phenomenon: 
a scalar field with a superluminal phase velocity considered in 
an appropriately Lorentz-boosted frame will have a vanishing 
kinetic term [20]. 
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FIG. 1. Conformal diagrams showing portions of de Sitter space charted by closed (left), flat (middle) and open (right) 
foliations. Thick lines indicate coordinate singularities. Superimposed are lines of constant isotropic time and radius for each 
foliation. 


metric describes a de Sitter space. We first discuss the 
conformal and three isotropic charts of the de Sitter space 
that play a prominent role in both the construction of 
solutions and the investigation of their global properties. 
We then turn to three families of self-accelerating solu¬ 
tions and discuss their global structure, namely the ap¬ 
pearance of their determinant singularities. 


A. De Sitter Charts 

One special feature of de Sitter space is that there is no 
preferred temporal coordinate to define a foliation with 
respect to. Sections of the full spacetime can therefore 
be charted by isotropic coordinates where the constant 
time slices have positive, negative or zero spatial curva¬ 
ture. The conformal diagram for de Sitter space can be 
constructed from the positive curvature (closed) foliation 
of the spacetime, where the line element takes the form 


the line element ( 22 ) takes its closed isotropic form 


ds 2 = -df 2 + 


cosh ( Ht c ) 


[l + (Hr c ) 2 /4\ 


(dr 2 + r 2 dfi 2 ) , (24) 


where t c £ (— 00 , 00 ), r c £ [0, 00 ). These coordinates 
chart the entire de Sitter spacetime. Similarly, defining 
the coordinates 


e Htf = - 


Hr f = 


cos X + cos ?7 
sin 77 

siny 


cos x + cos 77 
obtains the flat isotropic form 


(25) 


ds 2 = -At) + e 2Htf (dr) + r) dill) , (26) 

wher eHtf £ (— 00 , 00 ), Hrf £ [0, 00 ). These coordinates 
chart the upper left half of the conformal diagram 77 > 
X — 7T. Finally, the coordinate definition 


ds 2 = ( —\ (-dr/ 2 + dx 2 + sin 2 xdfi 2 ) , (22) 

\ H sin 77 ) ' 

with the dimensionless conformal time 77 £ (— 7 r, 0 ) and 
the comoving radial coordinate x G [0, 7 r]. Here H 2 = 
Agff/3. We use the ( 77 , x) conformal diagram throughout 
to represent the spacetime. 

Closed, flat, and open isotropic coordinates can al¬ 
ternately be used to foliate portions of de Sitter space 
and are useful in finding solutions to the background 
massive gravity equations and for investigating pertur¬ 
bations. With the transformations 

sinh (Ht c ) = — cot 77 , 

Hr c = 2tan(x/2), (23) 


In [tanh(7Tt 0 /2)l = tanh 1 (——, 

\cosxJ 

2tanh _ 1 (Ur 0 /2) = tanh -1 ( ———^ , (27) 

\ cos 77 ) 

gives the open isotropic form 


ds 2 = -df 2 


sinh ( Ht a ) 

Ll — (Hr 0 ) 2 /d\ 


(dr 2 


r 2 dfi 2 ) 


(28) 


where Ht a £ (0, 00 ), Hr a £ [0,2). These coordinates 
chart the upper left wedge of the conformal diagram 77 > 
X — 7t/2, corresponding to 1/8 of the space. We display 
these charts with the curves of constant tj and r t in Fig. 1. 

Self-accelerating background solutions appear superfi¬ 
cially different in different isotropic coordinates. More 
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form 

(/,* cot 2 % - / 2 ) sin 2 V + f 2 x + f, x f,v cot X sm(27?) 

= -f|csc V (30) 

Below we consider three explicit families of solutions 
to Eq. (30), corresponding to different self-accelerating 
backgrounds whose perturbations we will examine in the 
subsequent sections. 

As an aside, we note that (30) can be cast in a partic¬ 
ularly simple form, 


f2 _ f2 _ Hy 

At Ap ^2’ 

by a further coordinate transform [46] 

cosx 


P = 


sin g 


r = — cot g, 


(31) 


(32) 


FIG. 2. Constant unitary gauge radius g curves, given by 
Eq. (29). These are common to all bi-isotropic solutions. 


importantly, differences in the constant time surfaces 
cause the kinetic terms for the perturbations to appear 
differently in each case, as emphasized by Ref. [18]. Of 
course, the causal structure of the conformal diagram re¬ 
mains unchanged and serves to highlight the spacelike, 
timelike or lightlike nature of curves rather than coordi¬ 
nate dependent definitions of simultaneity. For this rea¬ 
son, it will be convenient to plot characteristics on the 
conformal diagram. 


B. Explicit Solutions 


We now combine the isotropic de Sitter charts with the 
formalism of §IIB to construct vacuum solutions to the 
massive gravity equations of motion. 

Radial solution: g and equations of motion 

Since closed isotropic coordinates chart the whole 
spacetime, it is convenient to use these coordinates to 
find self-accelerating solutions. Using the closed line el¬ 
ement Eq. (24) in the radial unitary gauge solution of 
Eq. (10), we obtain 


cosh (Ht c ) irnsiny 

9 = x o ,/tt J/A r <= = -' 

1 + [Hr c y /4 H sing 


(29) 


Fig. 2 shows the contours of constant g in the conformal 
diagram. As the conformal diagram makes obvious, g is 
4-fold symmetric in the de Sitter spacetime since g —> 
— 7 r — 77 and/or x —!► 7 r — % provide the same values. 

We then look for a choice of unitary time / that solves 
Eq. (17), which in the coordinates (22) takes the explicit 


with r £ (—oo, oo), p £ (—oo, oo). 

Open solution: f Q 

The simple ansatz that / is a function of g alone leads 
to the solutions first identified in Ref. [7] through an open 
slicing construction, where / is given by 

= cot ry. (33) 

In open slicing, the spacetime and fiducial metrics are 
bi-diagonal (because Eq. (18) is satisfied) and manifestly 
homogeneous and isotropic. This property of open slicing 
applies beyond the de Sitter solutions considered here to 
a general FRW spacetime, giving this solution a special 
status. As pointed out by Ref. [18], and rediscovered by 
Ref. [46], bi-diagonality does not hold in the closed or 
flat slicings of de Sitter. However, the ability to define 
multiple homogeneous and isotropic slicings and thread¬ 
ings of the spacetime is special to de Sitter. For a general 
FRW spacetime, this freedom no longer exists since the 
homogeneity of an evolving background density picks out 
a unique time slicing. 

Curves of constant unitary gauge time for this solu¬ 
tion are plotted in Fig. 3. Note that the coordinate pair 
( f 0 ,g ) maps to two different ( 77 , x) points in the space- 
time related by \ —^► tt — x and correspondingly the deter¬ 
minant of the Jacobian transformation to unitary gauge 
is zero along x = 7t/2. Therefore, more than one copy 
of the Minkowski fiducial space is required to cover the 
entire spacetime. For this solution, the singularity lies 
in the past light cone of an observer at x = 0 and hence 
boundary conditions that represent the continuation with 
a second fiducial metric are required. Note that—as men¬ 
tioned in § 11 — such a rule for continuation is ad hoc as it 
must be imposed by hand. 

Stationary solution: fc 

Self-accelerating solutions where the spacetime met¬ 
ric in unitary gauge is stationary were first identified in 
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FIG. 3. Constant unitary time curves for various well-known self-accelerating solutions: open f 0 (left), and stationary fc for 
C = 1/2 (middle) and fc for C = 1 (right). Solid lines indicate constant unitary time slices extending from x = 0 which end at 
a determinant singularity (red lines). Dashed lines indicate an extension of solutions with a second copy of the fiducial metric. 


Ref. [13]. In conformal coordinates, unitary time takes 
the form 


fc = 


Xp 

CH 


in 


C 2 e m f 

i -y 


V = ^/l + C 2 (sin 2 x/ sin 2 rj — 1 ), 


(34) 


where C £ (0,1] is a free parameter, y £ [0,oo), and 
tf(x, 7?) is the extension of flat isotropic time to the full de 
Sitter space using Eq. (25). The coordinate pair ( fc,g ) 
maps to two spacetime points ( 77 , x) and (— 7 r — 77 , 7t — x)- 
The inversion singularity thus is at x = and is char¬ 
acterized by a divergence in unitary time fc —t 00 which 
prevents its differentiability. Again, two copies of the 
Minkowski fiducial space are required to cover the space- 
time. In this case, the singularity is along the past light 
cone of the observer at x = 0 at the terminal confor¬ 
mal time 77 = 0. Conformal diagrams showing curves of 
constant fc for C = 1/2 and C = 1 are plotted in Fig. 3. 

A different construction of this solution that starts 
in static de Sitter coordinates [46] is in fact equivalent 
to Eq. (34) after relating their parameter q to C as 
q 2 = 1/C 2 — 1. The equivalence between the two classes 
of solutions can also be seen directly in the construc¬ 
tion itself since Ref. [13] showed that stationary unitary 
coordinates and static de Sitter coordinates are simply 
related by a radially dependent offset to the static time 
coordinate. 

New solution: f a b 

Finally, we can generalize a particular solution intro¬ 
duced in Ref. [18], to a new two parameter class, where 
the temporal Stiickelberg field takes the following form 


fab = (a - cot 77) 2 - {b - cos x esc rj) 2 . (35) 


Here, the parameters a, b can take on any real value 
except for a = b = 0. For that case, we have /q 0 = 
g 2 — ( Xq/H ) 2 , such that unitary time and radius cannot 
identify a unique spacetime point. Ref. [18] previously 
considered the case of a = 0 , b = 1. 

In the full spacetime f a b is not guaranteed to be real 
and so unitary gauge ends in a determinant singularity 
where W = ± 00 . Approaching this point, the change 
in unitary time (or proper time measured by the fiducial 
metric) per unit conformal time (or proper time mea¬ 
sured by the spacetime metric) diverges. Unlike the case 
of W = 0 , an analytic continuation of unitary coordi¬ 
nates would make a copy of the fiducial metric with a 
Riemannian rather than Lorentzian signature. Thus, un¬ 
like the other solutions, it does not appear that f a b can 
be continued beyond this type of determinant singular¬ 
ity with copies of the fiducial metric in the same class of 
solutions. 

Changing the sign of a reflects solutions vertically 
across 77 = — 7 t /2 whereas changing the sign of b reflects 
horizontally across x = 7 t/ 2 . The determinant singular¬ 
ities W = ±cc occur where f a b = 0 and bound regions 
past which the solutions cannot be continued within the 
class. These singularities intersect at 


cot 77 = a, cos x = — (36) 

V 1 + a z 

if \b\ < \J 1 + a 2 . The W = 0 singularity occurs along 
the curve a cos x = b cos 77 and intersects both W = ±00 
singularities at their crossing point. 

We show several examples from this class in Fig. 4. 
Displacing b from zero at a = 0 breaks the x or left-right 
symmetry of the conformal diagram allowing constant 
unitary time surfaces to foliate the spacetime near one 
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FIG. 4. Constant unitary time curves for the new f a b class of self-accelerating solutions. In addition to the zero determinant 
singularity (thick red lines) there also appear singularities where the determinant is infinite (blue lines), in the a = b = 1 case 
the two coincide. Across the infinite determinant singularity, a naive analytic continuation would require a Riemannian rather 
than Lorentzian second fiducial metric. 


pole with the other pole hidden behind a W = ±oo de¬ 
terminant singularity. The unbroken 77 or top-bottom 
symmetry has a corresponding W = 0 singularity at 
77 = — 7t/2 across which solutions can be extended with a 
second fiducial metric. Displacing a from zero at 6 = 0 
does the converse, allowing constant unitary time sur¬ 
faces to foliate the spacetime near both poles at either 
the top or bottom sections of the diagram with a W = 0 
singularity at x = vr/ 2 , again with a second branch of the 
solution across it. Displacing both equally makes one of 
the W = ±00 curves coincide with W = 0, allowing folia¬ 
tion around only one pole and only on the top or bottom 
section. 


C. Determinant Singularities 

All three families of solutions exhibit determinant sin¬ 
gularities. Worldlines of observers can intersect these 
singularities in the bulk of the spacetime, i.e., after only 
a finite amount of proper time has elapsed. For other ob¬ 
servers, the singularity lies in the past light cone. Both 
properties imply that the incompleteness of the solutions 
there cannot simply be ignored. 

In the open f a solution, a naive continuation past the 
W = 0 singularity keeps the Stuckelberg fields or uni¬ 
tary gauge coordinates continuous and differentiable but 
multivalued. In the static fc solution, naive continua¬ 
tion keeps them continuous but not differentiable at the 
singularity. In the new f a b solutions, naive continuation 
beyond the W = ±00 singularities is not possible and 
when they intersect at a specific point in the spacetime 
bulk they do so at the W = 0 singularity. 

Finally, by pushing the W = ±00 singularities of the 


fab case to the top or bottom of the conformal diagram it 
can be made to resemble the f 0 and fc solutions. The f 0 
solution picks out the same unitary time surfaces as the 
limiting case of a —> 00 , 6 = 0 since f Q and f a b are then 
linearly related. Likewise, the a = 6 —> 00 limit of f a b 
leads to the same surfaces of constant unitary time as the 
fc solution for C — > 0 in the upper right portion of the 
conformal diagram, although in this case f a b is a more 
general function of fc- The second copy of fc for C — > 0 
in the lower left corresponds similarly to a = b —» — 00 . 

As we shall see in the next section, the appearance of 
these singularities influences the characteristic curves on 
which information about perturbations propagate. 


IV. PERTURBATION CHARACTERISTICS 

In this section we use the method of characteristics to 
explicitly solve the equation of motion ( 20 ) for the field 
fluctuation 6 r in the three families of background solu¬ 
tion studied in §IIIB. The characteristic curves of this 
linear differential equation define hypersurfaces along 
which spherically symmetric stress-energy perturbations 
propagate. The conformal diagram of the curves exem¬ 
plifies the difference between superluminality, strong cou¬ 
pling, and the well-posedness of the Cauchy problem. 


A. Method of Characteristics 

We can solve the 6 T equation of motion for pertur¬ 
bations propagating on the background solutions using 
the method of characteristics. The characteristic curves 
of this equation also clarify the causal structure of solu- 
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tions. Transforming Eq. (20) from isotropic coordinates 
to the conformal coordinates ( 77 , %) leads to a differential 
equation of the form 


Vr,sr tV + v x dr >x + Adr = o. (37) 


Here, 


V v = f, x - tan x tan 77 /^, 

V x = f tT) + (cot % + esc 2 r) tan x) tan 77 f tX . (38) 

This equation depends on the particular background 
/(? 7 , x) around which we perturb and describes propa¬ 
gation of the <5T perturbations in the regions of de Sitter 
space where / is defined. While ^( 77 , x) is not important 
for the construction of characteristics themselves, it does 
enter into determining the field profile for dT along the 
characteristics. For completeness, it is given by 


N 


= sm ?7 


I V n p 
esc 77 \ ~ R 


where 


. 2 X 

+ COS - 


A* 


J ,v 

sec 2 'X ( w 


,2 X ( v x 


N 


J >x 


X 


R n = 2 cos x sin x cot - esc 77 , 

R x = n cos 2 ^ sin 2 x cot 77 esc 77 , 

„2 


N = esc x sec X t an V 


X 


2 fj, cos 77 tan ^ f x 


, 2 X ( . , x o 

+ sec - \jj, cos x sin 77 /,„ - — 


(39) 


(40) 


The most general solution to Eq. (37) can be found 
from its characteristics, which are integral curves for the 
vector field 


V = (V„V X ), (41) 

with the first component in the 77 direction. Information 
from boundary conditions on dT, specified on a surface 
that intersects the characteristic curves, propagates along 
those curves to provide the general solution in the bulk. 
Thus conformal diagrams of characteristic curves for par¬ 
ticular background solutions present a succinct way of 
describing their casual structure. In particular, the char¬ 
acteristic curves are 


timelike : |VJ,| > 

lightlike: \V v \ = \V x \, 
spacelike : |V^| < \V X \. (42) 

Since Eq. (37) is a first order differential equation, the 
characteristic analysis yields complete solutions for dT, 
including both smooth and discontinuous cases. This 
should be contrasted with a characteristic analysis of a 
second order differential equation where in most cases 
the characteristics only describe the propagation of dis¬ 
continuous fronts. Once a solution for <5T is specified from 
boundary data in the conformal coordinate system it can 
be expressed in any coordinate system since it transforms 
as a scalar function in spacetime. 


B. Explicit Solutions 

We now apply this formalism to the explicit back¬ 
ground solutions derived in §IIIB and identify the char¬ 
acteristic curves upon which isotropic stress-energy per¬ 
turbations propagate. 

Open solution: f Q 

For the f a background solution (33), the general solu¬ 
tion given by the characteristics is 

<5r o (77,x) = Sm 2 V F((j) 0 ), 
sin x 

MV,X) = cos x/sin 77 , (43) 

where F is a completely arbitrary function of its argu¬ 
ment (f> 0 . Curves of (j) 0 = const., shown in Fig. 5 (left), 
are the characteristics of the equation and comparison 
with Eq. (27) shows that they coincide with constant 
open time t a = const, surfaces in the wedge charted by 
the open coordinates. Correspondingly, in open isotropic 
coordinates, Eq. (20) loses its kinetic term (A to = 0) in¬ 
dicating that initial conditions at t 0 = const, do not prop¬ 
agate off that surface. Open time slices therefore are dis¬ 
connected, leading to an unspecified evolution of <5T. In 
particular, information propagates instantaneously along 
the characteristics as measured by open time. This be¬ 
havior does not violate local conservation of energy and 
momentum since they are both zero for any dT 0 in these 
coordinates [17]. 

In flat or closed coordinates, these features take a su¬ 
perficially different form. Characteristics are still space¬ 
like but the finite stress components in open coordi¬ 
nates transform to energy and momentum. This energy- 
momentum is conserved by the solution (43) as a conse¬ 
quence of Eq. (20). As discussed in §IIC, this conserva¬ 
tion law is the transformation of the constraint on the 
spatial stresses in open coordinates. 

On the other hand, the stress-energy emanates from 
X = 0 where the spatial profiles of the dr o solutions 
diverge. Although it is not manifest in the linearized 
treatment from the quadratic Lagrangian, generically one 
might expect perturbation theory to break down here 
with higher order interactions making the fluctuations 
strongly coupled. This class of solution was excluded 
by the analysis of Ref. [47] by implicitly demanding reg¬ 
ularity at x = 0 , leading to their conclusion that the 
Stiickelberg fields contained no degrees of freedom. Our 
characteristic analysis makes this assumption explicit: 
it corresponds to assigning boundary conditions on the 
timelike x = 0 curve. More generally, there is no space¬ 
like Cauchy surface at all in the open wedge, since a 
surface must be timelike to intersect all characteristics. 
The lack of a spacelike Cauchy surface is diffeomorphism 
invariant. Thus, even in flat or closed slicing, counting 
degrees of freedom according to initial data would lead to 
the conclusion that the f 0 case does not possess a degree 
of freedom that admits a well-posed Cauchy problem. 
The spacetime defined by the solution f a is an example 
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FIG. 5. Characteristic curves of <5r for the background solutions shown in Fig. 3. For <5F 0 , characteristics coincide with constant 
open time curves in the open wedge of Fig. 1 and no spacelike surface intersects all characteristics. For dFj/ 2 , characteristics 
are also spacelike but those within the past light cone at \ = 0 do all intersect the spacelike surface y = — n. For the special 
case of <5Fi, the characteristics are all lightlike. Red lines divide the characteristics on either side of the determinant singularity. 


of a spacetime that is not globally hyperbolic . 2 

While imposing a regular boundary at y = 0 for these 
characteristics might seem reasonable in the open wedge 
itself, it is interesting to track characteristics intersect¬ 
ing this boundary in different parts of the de Sitter space. 
Fig. 5 (left) shows that in the lower left wedge, where the 
characteristics are mirror images of the open wedge, there 
is a spacelike Cauchy surface at y —> —i r despite superlu- 
minality, but the characteristics end rather than begin at 
X = 0. Stress-energy on a characteristic then flows into 
the origin requiring a boundary condition or nonlinear 
completion of the theory to specify its further evolution. 
Thus non-singular behavior on the open wedge at y = 0 
may require special initial conditions in the lower left 
wedge. On the other hand this behavior is analogous to 
collapse of a perfectly spherically symmetric density shell 
in ordinary linearized, Newtonian gravity. The finite an¬ 
gular momentum of generic perturbations may prevent 
singular behavior—we leave investigation of this possi¬ 
bility to future work. In any case, while the inwardly 
directed characteristics might signal strong-coupling at 
y = 0 where the fluctuations diverge, like a black-hole 
nothing emanates classically from this point within the 
lower left wedge. 

Interestingly, in the central diamond the <5T 0 character¬ 
istics are timelike, or subluminal, and admit a well-posed 
Cauchy problem. However, the characteristics are split 
by the determinant singularity into left and right halves. 


2 The existence of such a pathological spacetime as a solution to 
dRGT does not imply pathologies inherent to the theory itself; 
indeed, general relativity admits solutions which are not globally 
hyperbolic. 


Although information in <5T 0 never crosses this singular¬ 
ity in the linearized theory, the worldlines of other par¬ 
ticles can. From the perspective of a massive gravity 
theory with only one copy of the fiducial metric, the de¬ 
terminant singularity appears as another spatial bound¬ 
ary condition. Even allowing the theory to be extended, 
the boundary condition must be imposed by hand at each 
order in perturbation theory, e.g ., by demanding that the 
Stiickelberg fields are smooth across the boundary to the 
second copy of the fiducial metric [43]. 

Stationary solution: fc 

The characteristics of perturbations around the fc 
backgrounds share some, but not all, of the properties 
of the f 0 background. Here Eq. (37) is solved by 


sr c(v,x) 


sin 2 r) F{<j) C ) 
sin 2 x y(r], x) ’ 
sinr?[l - yfa y)] 
cos x + cos 77 


(44) 


where F is again an arbitrary function of its argument 
and y was defined in Eq. (34). In Fig. 5 (mid), we show 
the 1 fic = const, curves for C = 1/2. These charac¬ 
teristics are spacelike across the whole spacetime. Even 
though (f>c = const, does not coincide with the open, flat 
or closed choice of isotropic time, these surfaces foliate 
the spacetime. Thus in principle we could define a new 
choice of time tc = 4>c for which the kinetic term for 
<5Tc vanishes, just like it does for <5T 0 in the open slicing. 
With this choice of time coordinate, information along 
the characteristics again propagates instantaneously and 
energy-momentum conservation becomes a spatial con¬ 
straint. 

On the lower left half of the 5Tc diagram, the causal 
structure resembles the lower left wedge of the (7T 0 case. 
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FIG. 6. Characteristic curves of <5r a b for the new class of f a b solutions illustrated in Fig. 4. The determinant singularities 
strongly influence the causal structure of the characteristics. These backgrounds share many of the features of the f 0 and fc 
spacetimes discussed in Fig. 5 and provide new ones due to the W = ±oo determinant singularities (blue curves). 


Here ij —ir provides a spacelike Cauchy surface and the 
characteristics end at y = 0 with a divergent ST 1/2 radial 
profile. On the upper right half, the structure resembles 
the open wedge of <5r o with characteristics emanating 
from the other pole of the closed de Sitter space. How¬ 
ever, the two halves are divided by the null line that des¬ 
ignates the determinant singularity where / is not contin¬ 
uously differentiable. Thus the outgoing characteristics 
of upper right half are not within the past light cone at 
X = 0 unlike the dT 0 case. This example illustrates that 
superluminality, a well-defined Cauchy problem, and sin¬ 
gular field profiles in the past that might indicate prob¬ 
lems originating from a point of strong coupling are not 
one-to-one related. 

In the special case of C = 1, y = — sin y/ sin rj and 
the characteristics of <5T 1 become null (see Fig. 5 and 
Ref. [17]). This special case was interpreted by Ref. [13] 
as having no vector in the background, based on a de¬ 
coupling limit analysis. Although this choice eliminates 
the superluminality found in other solutions, the lightlike 
determinant singularity still exists for this case. While 
it only occurs on the past liglitcone of the observer at 
X = 0, 77 = 0, it would be within the past light cone 
at other points in the spacetime and still require null 
boundary conditions to define global solutions. Likewise 
other observers in the lower left region would reach the 
determinant singularity after finite proper time. Even 
with solutions continued beyond the singularity these ob¬ 
servers would begin to cross characteristics for which no 
spacelike Cauchy surface exists. 


lution is 


5Fab(V,X) 

0ab(V,X) 


esc 2 x si 1 ! 3 V 
cos r] — a sin rj 
cos x — b sin 77 
cos rj — a sin r\ 


F(<t>ab), 


(45) 


This solution serves to illustrate that the nature of the 
characteristics is strongly influenced by the global struc¬ 
ture of the determinant singularities discussed in §IIIB. 
Characteristics run tangent to these singularities as in¬ 
formation in the perturbations cannot cross these curves 
(see Fig. 6 ). If | 6 | < \/l + a 2 , the W = ±oo and W = 0 
singularities intersect within the spacetime bulk at a 
point given by Eq. (36). Though the characteristics tan¬ 
gent to these singularities also intersect, the perturbation 
amplitude <5r a b diverges at this point and hence marks a 
point where strong coupling is likely. 

Null lines that emanate from this point serve as sep- 
aratrices which divide regions of super- and subluminal 
propagation much like the f a case. This family of exam¬ 
ples shows that these null lines are not necessarily related 
to special curves in the open, closed or flat de Sitter slic¬ 
ing. Where these null lines intersect the poles, the char¬ 
acteristics change from ingoing into the pole, which can 
be all be intersected by a single spacelike Cauchy surface 
in the past, to outgoing which cannot. For the special 
case that |a| = | 6 | one of these null lines coincides with 
the determinant singularity. 


C. Superluminality and Bi-Isotropy 


New solution: f a b 

For the new class of f a b backgrounds, the general so- 


In all but one example (<5rc with C = 1), the charac¬ 
teristics intersect the poles y = 0 , 7 r at right angles at all 
but a finite number of points. In conformal coordinates, 
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perturbations propagate instantaneously there. As we 
now show, this is a generic property of bi-isotropic solu¬ 
tions. In our explanation we focus only on the % = 0 pole 
as treatment of the other pole is analogous. Expanding 
unitary time around x = 0 i n a series 

f(V, X) = h 0 (r]) + hi(rj)x + h 2 (y)x 2 + 0{x 3 ) (46) 

and using the x expansion of the equation of motion (30), 
one can show that hi is always identically zero. Physi¬ 
cally this means that the bi-isotropy condition and the 
ability to make arbitrary Lorentz boosts in the fiducial 
space ensures that unitary time and conformal time can 
be locally aligned. The difference between surfaces of uni¬ 
tary time and conformal time thus grow at most quadrat- 
ically with the distance to the pole x 2 or more generally 
differ only by the curvature corrections to a locally flat 
spacetime metric. Since the same is true between confor¬ 
mal time and closed, flat, or open isotropic time, instan¬ 
taneous propagation in conformal coordinates implies in¬ 
stantaneous propagation in isotropic coordinates at the 
poles. In fact, the arguments below for the generality of 
spacelike characteristics apply beyond the vacuum cases 
considered here to spacetime metrics that are locally flat 
near the bi-isotropy point. 

Expanding the vectors defined by Eq. (38) in the same 
fashion, 


( 0 { X ) \ 

\ 2h 2 tan r] +h 0 ,r,+0(x) )' 


(47) 


In a neighborhood of a general point ( 770 , 0) the leading 
order terms are finite and the vector field V is aligned 
with the curves of constant conformal time 77 , signaling 
instantaneous propagation. Using the expansion of the 
background equation of motion (30), it is possible to re¬ 
late h 2 and ho,n to prove that infinite superluminality at 
the given point is avoided if and only if 


ho, v (v 0 ) = ± 


x 0 

H sin 770 


(48) 


As seen with our examples, this condition can be satisfied 
for all 77 (<5Tc with C = 1) or at a discrete set of points 
(<5r o and typical <5r a b) . 

In fact, for the special points 770 at which this condition 
is satisfied, the characteristics are always luminal, never 
subluminal. Expanding / to third order in x allows V 
to be expanded into linear order in both x an( 4 V ~ Vo- 
With repeated use of the equation of motion (30) we then 
find that each characteristic curve hitting x = 0 at 770 is 
luminal. For h 3 (r]o) = 0 there are two such characteris¬ 
tic curves—one incoming and one outgoing—forming the 
typical luminal separatrices we see in Fig. 5 (left) and 6 . 
Even for these cases, the characteristics are superlumi¬ 
nal at all but a finite number of points. For h 3 (rjo) 7 ^ 0, 
there is only one characteristic which is either incoming 
or outgoing, depending on the sign of ho, v /h 3 . The (7=1 
example of Fig. 5 (right) exhibits a unitary time solution 


/1 with these very special properties. Together with so¬ 
lutions trivially related by / —► —f and/or 77 —» —n — 77 , 
it is in fact the unique bi-isotropic vacuum solution that 
evades superluminality entirely as can be shown by inte¬ 
grating the equation of motion (30) from the boundary 
conditions (48) along null coordinates. 


V. DISCUSSION 

We have investigated the causal structure of three 
families of vacuum solutions in dRGT massive gravity. 
In particular, we have constructed the conformal dia¬ 
gram of characteristic hypersurfaces for isotropic stress- 
energy perturbations around these backgrounds by ex¬ 
ploiting the first-order structure of their equation of mo¬ 
tion. These examples provide fertile ground for study¬ 
ing the interplay between superluminality, an ill-posed 
Cauchy problem, and strong coupling as well as issues 
that arise for the global structure of spacetime in a the¬ 
ory with two metrics. 

The f 0 solution of Ref. [7] manifests aspects of all of 
these issues. This solution is distinguished because in 
open slicing both the background spacetime and fiducial 
metric are simultaneously homogeneous and isotropic. In 
this slicing, the kinetic term of perturbations vanishes, 
which was taken as an indication of possible problems 
with strong coupling [47]; this is supported by instabili¬ 
ties identified around anisotropic backgrounds [48], and 
the lack of an isotropic degree of freedom identified by 
a Hamiltonian analysis [18]. However, its kinetic term 
and Hamiltonian degrees of freedom only vanish in open 
slicing [18]. 

Our characteristic analysis clarifies these issues. Ab¬ 
sence of a kinetic term in the quadratic action in a par¬ 
ticular choice of slicing indicates instantaneous propaga¬ 
tion of perturbations in the given frame and more gener¬ 
ally perturbative superluminality or spacelike character¬ 
istics in any frame. The Hamiltonian analysis fails in the 
pathological choice of slicing along characteristics since 
the Hamiltonian is no longer associated with time evolu¬ 
tion. Energy-momentum conservation likewise becomes 
a spatial constraint equation, confusing the counting of 
degrees of freedom. While it is the easiest of the inter¬ 
related issues to diagnose, it does not necessarily indi¬ 
cate strong coupling nor an ill-posed Cauchy problem. 
Through a particularly poor choice of coordinates, it is 
possible to make the kinetic term for superluminal fluc¬ 
tuations disappear even in a free theory. Likewise, in¬ 
formation cannot propagate forward from a pathological 
choice of initial value surface. 

On the other hand, for the f 0 solution, the characteris¬ 
tic analysis does additionally reveal an ill-posed Cauchy 
problem. Since spacelike characteristics originate from 
the spatial origin in the open chart of de Sitter, there 
is no spacelike surface that intersects all characteris¬ 
tics. Hence this particular spacetime does not admit a 
well-defined initial value problem, as it lacks a spacelike 
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Cauchy surface. In contrast to superluminality, nonex¬ 
istence of spacelike Cauchy surface always implies the 
quadratic action is pathological. Further, in contrast to 
the Hamiltonian identification of degrees of freedom, this 
concept and its relation to identifying degrees of freedom 
by the amount of initial data required is diffeomorphism 
invariant. 

Finally, for the f a solution, finite perturbations speci¬ 
fied along characteristics diverge in amplitude at the spa¬ 
tial origin from which they emanate in the open chart. 
Given that the theory contains nonlinear interactions, 
this indicates that perturbation theory almost certainly 
breaks down leaving the theory strongly coupled there. 
Hence in this case, the ill-posed Cauchy problem likely 
originates from a point of strong coupling. More pre¬ 
cisely, strong coupling occurs when the effective field 
theory breaks down due to ever higher order interaction 
terms becoming important. While strong coupling can¬ 
not be diagnosed from the quadratic action alone, the 
divergence or discontinuity of perturbative solutions at a 
given spacetime point is a good indicator of strong cou¬ 
pling there. 

Vanishing of kinetic terms is often used as a related 
proxy for strong coupling: once variables are canonically 
normalized, interaction terms pick up negative powers of 
the coefficient in front of the kinetic term, which drives 
the effective strong coupling energy scale to zero when the 
kinetic term vanishes. However, our example shows that 
without examining the higher order terms themselves, 
one cannot immediately determine whether the vanishing 
of kinetic terms instead simply arises from a poor choice 
of coordinates. 

The fc, or stationary class of solutions, serves to fur¬ 
ther distinguish these concepts. Here generic solutions 
also admit spacelikc characteristic curves and hence su- 
perluminality. By choosing a time coordinate that is or¬ 
thogonal to these spacelike surfaces, we again recover a 
representation where kinetic terms and dynamics vanish 
in favor of spatial constraints. Unlike the open f a solu¬ 
tion, all characteristics of this class that are within the 
past light cone of an observer at the spatial origin in¬ 
tersect a spacelike surface. Hence the Cauchy problem 
is well-posed for this family of solutions and observers. 
Furthermore, the characteristics end rather than begin at 
the spatial origin where the field perturbations diverge. 
While strong coupling at the origin is again likely, this 
distinction is crucial. As in spherical collapse of pertur¬ 
bations in general relativity (with Newtonian gravity as 
the corresponding effective theory), the formation of a 
singularity at the origin does not necessarily invalidate 
the effective theory far from the singularity if nothing 
escapes it classically. 

Finally, the new class of f a b solutions serves to high¬ 
light issues with the global structure of the spacetime 
that can occur in theories with two metrics. In all three 
families of solutions, the fiducial Minkowski space does 
not cover the whole de Sitter spacetime. Specifically, 
there are points where the diffeomorphism-invariant ratio 


of the determinants of the physical and fiducial metrics 
vanishes, diverges or is undefined. Here the spacetime 
becomes geodesically incomplete. While the f a and fc 
solutions allow straightforward but ad hoc analytic con¬ 
tinuation past these determinant singularities to define 
a solution in the whole spacetime, the f a & ones do not. 
Furthermore, since characteristics do not cross these sin¬ 
gularities, continuity must be imposed not just for the 
background but also for the perturbations. 

The superluminality exhibited in all three classes of 
solutions is a necessary condition of the bi-isotropic con¬ 
struction except in one unique case. Generically, starting 
at the bi-isotropy point characteristics are superluminal 
across separations comparable to the spacetime curva¬ 
ture as a consequence of the alignment between the lo¬ 
cally flat spacetime metric with the fiducial flat metric. 
Exceptional cases produce luminal but never subluminal 
characteristics that intersect this worldine. For the sin¬ 
gle case of the stationary fi solution, all characteristics 
are strictly luminal. Unlike the usual second-order sys¬ 
tem where characteristics define only a front velocity of 
discontinuous solutions, the characteristic analysis of our 
first-order system is fully general. In particular, it al¬ 
lows the construction of smooth wavepackets that also 
propagate superluminally and do not on their own imply 
strong-coupling. 

While these examples do serve to distinguish the con¬ 
cepts of superluminality, spacelike Cauchy surfaces and 
strong coupling, it is important not to overinterpret 
their consequences for the dRGT theory itself. They 
represent examples of tree-level or classical propagation 
of perturbations on specific, perhaps pathological, self- 
accelerating backgrounds. Superluminality in the full 
theory would be in tension with the dRGT theory ad¬ 
mitting a local and Lorentz-invariant UV completion [20] . 
We emphasize, however, that the characteristic analysis 
alone is insufficient to establish whether or not superlu¬ 
minal propagation is truly present—in order to discern 
this, it has been claimed that one must know about the 
Green’s function at infinite frequency [21-23], which goes 
beyond the classical approximation that the characteris¬ 
tic analysis entails. 

It is also worth mentioning that the presence of super¬ 
luminality does not in and of itself imply acausality [41]. 
Indeed, in order to establish acausality, one would need 
to use the superluminal propagation to construct a closed 
timelike curve (e.<?., by using a point at which character¬ 
istic curves cross). Interestingly, we find (at least in the 
highly symmetric case we consider) that characteristic 
curves only intersect at singularities, indicating that such 
perturbations remain causal within the effective theory, 
even if we equate perturbative superluminality with true 
superluminality. It has been conjectured in Ref. [42] that 
this might be true in general: the chronology protection 
conjecture asserts that traversing a closed timelike curve 
cannot be achieved within the regime of validity of the 
effective theory. 

Finally in all of our examples, the Hamiltonian asso- 
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ciated with isotropic perturbations is unbounded [18]. 
While it is not clear how to interpret the Hamiltonian 
in these cases where the equations of motion are first or¬ 
der and there exists pathological choices of slicing where 
it is no longer associated with time-evolution, this may 
indicate that the features uncovered here are associated 
with an unstable background rather than generic conse¬ 
quences of the dRGT theory. 
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